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ABSTRACT 

We develop atmosphere models of two of the three Kepler-field planets that were known 
prior to the start of the K epler mission (HAT-P-7b and TrES-2). We find that published Ke- 
pler and Spitzer data for HAT-P-7b appear to require an extremely hot upper atmosphere on 
the dayside, with a strong thermal inversion and little day-night redistribution. The Spitzer 
data for TrES-2 suggest a mild thermal inversion with moderate day-night redistribution. 
We examine the effect of nonequilibrium chemistry on TrES-2 model atmospheres and find 
that methane levels must be adjusted by extreme amounts in order to cause even mild 
changes in atmospheric structure and emergent spectra. Our best-fit models to the Spitzer 
data for TrES-2 lead us to predict a low secondary eclipse planet-star flux ratio ( ^ 2x f 0~ 5 ) 
in the Kepler bandpass, which is consistent with what very recent observations have found. 
Finally, we consider how the Kepler-hand optical flux from a hot exoplanet depends on the 
strength of a possible extra optical absorber in the upper atmosphere. We find that the 
optical flux is not monotonic in optical opacity, and the non-monotonicity is greater for 
brighter, hotter stars. 

Subject headings: planetary systems - radiative transfer - stars: individual HAT-P-7, TrES-2 



1. INTRODUCTION 

Extrasolar planets are being discovered at an 
increasingly rapid pace: roughly a quarter of the 
currently known exoplanets (numbering more than 
450, as of June, 20f 0) were found since the begin- 
ning of 2009. 1 Still, only ^80 of the known planets 
have been seen to transit their parent stars. Tran- 
sits break the degeneracy between mass and incli- 
nation, they allow direct measurement of planetary 
radii, and they make possible precise measurements 
of planetary fluxes from secondary eclipse observa- 
tions. The Kepler mission, which is predicted to 
find many new transiting planets, is, therefore, par- 
ticularly exciting. 

Three transiting planets in the Kepler field were 
identified prior to the beginning of the Kepler mis- 
sion - TrES-2, HAT-P-7b, and HAT-P-lfb (occa- 
sionally referred to as Kepler-fb, -2b, and -3b). 
Spiegel et al. (2010) have already published a range 
of possible atmospheric models of HAT-P-llb; here, 
we consider HAT-P-7b and and TrES-2. 

The InfraRed Array Camera (IRAC) instrument 
on the Spitzer Space Telescope has been a boon 
to exoplanetary science, providing observations that 
are diagnostic of atmospheric temperature and corn- 
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position for more than a dozen planets. It had four 
photometric channels, centered at 3.6 /im, 4.5 ^tm, 
5.8 /xm, and 8.0 /im. Recently, Christiansen et al. 
(2010) and O'Donovan et al. (2010) used IRAC to 
measure infrared fluxes from HAT-P-7b and TrES- 
2, respectively. 

HAT-P-7b, discovered by Pal et al. (2008), orbits 
a large, hot star (1.84i?©, spectral type F8). It is 
one of the most highly irradiated known explanets, 
with a substellar flux of ~4.8 x 10 9 erg cm~ 2 s -1 . 
Its orbit is significantly misaligned from the stellar 
spin vector, indicating a possible third body in the 
system (Winn et al. 2009; Narita et al. 2009). It 
is a particularly interesting object in part because 
Kepler's exquisite photometry has allowed measure- 
ment of ellipsoidal variations in the star induced by 
the planet's tidal field (Welsh et al. 2010). TrES- 
2, by contrast, orbits a nearly solar-type star in a 
nearly grazing orbit (O'Donovan et al. 2006; Hol- 
man et al. 2007; Sozzetti et al. 2007; Raetz et al. 
2009). Mislis & Schmitt (2009) find a reduction in 
transit duration of ^3 minutes since 2006, and at- 
tribute this shortening to a change in inclination, 
although analysis by Rabus et al. (2009) does not 
corroborate such a large change. 

The atmosphere modeling strategy that we em- 
ploy here differs from several others that have been 
used. Similar to Barman et al. (2005) and Fort- 
ney et al. (2006), we calculate radiative equilib- 
rium, chemical equlibrium models. In contrast, 
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both Madhusudhan & Seager (2009) and Tinetti 
ct al. (2005, 2007) adjust both chemistry and ther- 
mal structure in order to find a best fit to the avail- 
able data, eschewing equilibrium solutions. This 
latter method produces chemical and thermal pro- 
files that are not the result of ab initio calcula- 
tions, but that might reveal non-equilibrium behav- 
ior. Yet another approach is taken by Showman 
et al. (2009) and Burrows et al. (2010), who sim- 
ulate three dimensional structure and dynamics in 
planetary atmospheres; so far this more sophisti- 
cated approach has not produced better fits to ob- 
servations than the one-dimensional radiative mod- 
els described above. 

Hubeny et al. (2003) were the first to suggest 
that an extra optical absorber in a hot exoplanet's 
upper atmosphere could lead to a thermal inver- 
sion. Observations seem to suggest such inversions 
(Hubeny et al. 2003; Fortney et al. 2006, 2008; Bur- 
rows et al. 2007, 2008a; Richardson et al. 2007; 
Spiegel et al. 2009; Madhusudhan & Seager 2009; 
Knutson et al. 2008, 2010), and the strong opti- 
cal absorber titanium oxide (TiO) has frequently 
been suggested as the possible culprit responsible 
for the inferred inversions. However, in the absence 
of strong mixing processes, the molecular weight 
of TiO would tend to make it settle to the bot- 
tom of the atmosphere. Furthermore, cold traps 
on the nightside and below the hot upper atmo- 
sphere on the dayside can cause TiO to condense 
into solid grains, which necessitates even stronger 
macroscopic mixing to keep TiO aloft in the radia- 
tively important upper atmosphere. Since the pho- 
tospheres are above the radiative-convective bound- 
aries, they are stably stratified; it is not obvious 
whether such strong mixing obtains in such a stable 
region (Spiegel et al. 2009). Various authors have 
tried to estimate the amount of macroscopic mixing 
and have found that, in some regions of the day-side 
atmosphere, the mixing might be vigorous enough 
to maintain TiO at altitude (Showman ct al. 2009; 
Li & Goodman 2010; Madhusudhan & Seager 2010), 
though Youdin & Mitchell (2010) point out that tur- 
bulent diffusivity in excess of 10 7 cm 2 s _1 might lead 
to overinflation of some planets through downward 
transport of entropy. Zahnle et al. (2009) suggest 
that sulfur photochemistry provides another avenue 
for achieving the additional upper-atmosphere op- 
tical opacity that is needed to produce hot upper 
atmospheres and inversions. Knutson et al. (2010) 
find that the presence of thermal inversions appears 
to be inversely related to the host stars' ultraviolet 
(UV) activity, raising the possibility that strong in- 
cident UV destroys molecules that may be responsi- 
ble for thermal inversions. Here, we remain agnostic 
on the matter and include, for modeling purposes, 
an extra source of optical opacity of as-yet unknown 
origin. 

This paper is structured as follows. In §2, we 
describe our ID atmosphere modeling strategy. In 



§3 and §4, we present our models of HAT-P-7b and 
TrES-2. We point out the effects that scattering and 
nonequilibrium chemistry could have on our models. 
In §5, we discuss the non-monotonic relationship be- 
tween the optical opacity in a planet's upper atmo- 
sphere and the planet's optical emergent radiation. 
Finally, in §6, we conclude, and in an Appendix, we 
discuss the various methods used in the literature to 
represent day-night redistribution in this modeling 
context. 

2. ATMOSPHERE MODELING 

As in our other recent studies, we use the code 
C00LTLUSTY (Hubeny et al. 2003; Sudarsky et al. 
2003; Burrows et al. 2006, 2008a; Spiegel et al. 2009, 
2010), a variant of the code TLUSTY (Hubeny 1988; 
Hubeny & Lanz 1995), to calculate radiative equi- 
librium irradiated atmosphere models. Our atomic 
and molecular opacities are generally calculated as- 
suming chemical equilibrium with solar elemental 
abundances (Sharp & Burrows 2007; Burrows & 
Sharp 1999; Burrows et al. 2001, 2002, 2005), al- 
though we also calculate a few nonequilibrium mod- 
els, described in §4. The irradiating spectra in our 
models are taken from Kurucz (1979, 1994, 2005), 
interpolated to the effective temperatures and sur- 
face gravities of HAT-P-7 and the host star of TrES- 
2. 

In addition to the observationally measured pa- 
rameters of our models (orbital semimajor axis, 
planet and stellar radii, planet and stellar surface 
gravities, stellar effective temperature), and in ad- 
dition to calculating equilibrium chemistry and ra- 
diative transfer, several other physical processes go 
into our models. These include Rayleigh scatter- 
ing (Sudarsky et al. 2000; Lopez-Morales & Seager 
2007; Burrows et al. 2008b; Rowe et al. 2008), heat 
redistribution, and the possible presence of an extra 
optical absorber that could explain the hot upper 
atmospheres and thermal inversions that have been 
inferred from infrared observations of several tran- 
siting planets. In particular, there are two key free 
parameters that we vary: P n and n' . 

• P n quantifies the efficiency of day-to-night 
heat redistribution, and is equal to the frac- 
tion of incident day-side heating that is reradi- 
ated from the nightside (Burrows et al. 2006, 
2008a). In the models in this paper, the re- 
distribution takes place between 0.01 and 0.1 
bars. 

• k' is an ad hoc extra source of optical ab- 
sorption opacity in the upper atmosphere, and 
is motivated by the thermal inversions that 
have been inferred from the infrared spectra 
of many exoplanets. n' is similar to the n e of 
some of our recent work (e.g., Burrows et al. 
2008a), except that, rather than a gray op- 
tical absorber, it has the same parabolic de- 
pendence on frequency as the corresponding 
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extra absorber in Lopez-Morales et al. (2009) 
and Burrows et al. (2010). 

The models whose properties arc plotted in Figs. 1- 
3 are summarized in Table 1. In addition to these, 
we also calculate several models with a modification 
to the code that allows an ad hoc extra source of 
optical scattering opacity, analogous to «/, but for 
scattering instead of absorption. 

3. HAT-P-7b 

HAT-P-7b is a 1.78M/, lMRj planet in an 
approximately circular orbit 0.0377 AU from the 
1.47 Mq, 1.84 R Q , F6 star HAT-P-7 (Pal et al. 
2008). In addition to being in the Kepler field, HAT- 
P-7b is interesting because Rossiter-McLaughlin 
measurements indicate that it is in a polar or retro- 
grade orbit (Winn et al. 2009; Narita et al. 2009). 
More relevant to modeling its atmosphere is the fact 
that, due to its close proximity to a large, relatively 
hot (^6350 K) star, HAT-P-7b experiences unusu- 
ally high stellar irradiation (<~4.8 x 10 9 erg cm~ 2 s _1 
at the substellar point). 

As an early confirmation that Kepler was per- 
forming well, Borucki et al. (2009) published 10 days 
of commisioning-phase data on HAT-P-7b. These 
data reveal a surprisingly large secondary eclipse 
depth in the Kepler band (^0.43-0.83 /jm), as the 
corresponding planet-star flux ratio decreases by 
(1.3 ± 0.1) x 10~ 4 when the planet passes behind 
the star. 2 Borucki et al. (2009) suggest that such a 
large contrast ratio could imply that the atmosphere 
absorbs strongly and has minimal redistribution to 
the night side (low P n , in our language). They es- 
timate a day-side temperature of 2650 ± 100 K. 

Infrared data for HAT-P-7b became available 
shortly thereafter, when Christiansen et al. (2010) 
presented secondary eclipse observations of the 
planet employing the IRAC instrument on the 
Spitzer Space Telescope. They found planet-star 
flux-ratios of (9.8±1.7) x 10~ 4 , (15.9±0.22) x 10~ 4 , 
(24.5 ± 3.1) x 10~ 4 , and (22.5 ± 5.2) x 10~ 4 at the 
IRAC 3.6-/im, 4.5-/im, 5.8-^m, and 8.0-^m, chan- 
nels, respectively. 

Although, as Christiansen et al. (2010) ar- 
gue, HAT-P-7b's atmosphere is irradiated strongly 
enough that it is probably too hot for the con- 
densates that might otherwise be expected to con- 
tribute significant scattering opacity, we neverthe- 
less tried adding an ad hoc extra scatterer to the 
upper atmosphere of some models. We find that 
the optical point can easily be matched by a model 
with significant optical scattering, but such models 
drastically underpredict the Spitzer data. Since the 
Borucki et al. (2009) speculation of an extremely hot 

2 This contrast ratio is similar to the expected dip in flux 
for an Earth-sized planet passing in front of a Sun-like star, 
and therefore indicated that Kepler should be capable of per- 
forming the mission for which it was designed. 



upper atmosphere with low albedo came before the 
infrared data were published, it might have been a 
little bit premature. On the other hand, the expec- 
tation of low albedo is well motivated by the high 
stellar irradiation, and seems to be confirmed by the 
infrared observations. 

In Table 1, we present five thermochemical equi- 
librium models of HAT-P-7b, four having extremely 
hot upper atmospheres (ranging from H2's ^3040 K 
to H5's ^3180 K), and a comparison model without 
a thermal inversion (HI). These five models span 
a (small) range of values of both P n (the degree 
of day-night redistribution) and k' (the strength of 
an ad hoc extra absorber, where the values are in 
cm 2 g _1 ). Models with values of P n larger than 0.1 
are not displayed in this table because such models 
underpredict the optical data. 

The five HAT-P-7b models of Table 1 are dis- 
played in Fig. 1. The top- left and top-right panels 
portray the wavelength-dependent planet-star flux 
ratios in the optical and the infrared, respectively, 
and the corresponding data (Borucki et al. 2009; 
Christiansen et al. 2010) are superimposed. In both 
panels, planet and star fluxes are integrated over the 
relevant bandpasses, and the resulting integrated 
planet-star ratios are displayed as solid filled circles. 
The bottom panel shows temperature-pressure pro- 
files for these five models. Models H2-H5, with ex- 
tra optical opacity, have strong thermal inversions 
in which the upper atmosphere is heated to tem- 
peratures <; 1500 K greater than they would be in 
the absence of the extra absorber (as represented 
by model HI). 3 The near- and mid-infrared spectra 
for inverted models H2-H5 are nearly indistinguish- 
able from one another, and are all reasonable fits 
to the IRAC data (although the models are ^l-2er 
higher than the data at 3.6 /im and 4.5 /an). Taking 
into account the optical data (top-left panel) , model 
H5 is the best fit to the available data. HAT-P-7b's 
real atmosphere is probably not represented by a ID 
radiative equilibrium model, such as model H5, but 
the available data suggest that significant flux might 
be absorbed high in the atmosphere, where the ra- 
diative timescale is short (Iro et al. 2005; Showman 
et al. 2008), which would imply that day-side heat 
would probably be reradiated before much advective 
redistribution has occurred. The association of high 
extra opacity with low day-night redistribution in 
the best-fitting model, therefore, is consistent with 
what should be expected. 

Christiansen et al. (2010) analyzed the thermo- 
chemical implications of the Kepler/ Spitzer data, as 
well. Using the method described in Madhusudhan 
& Seager (2009) and the four IRAC data points 
in conjunction with the Kepler data, they find 

3 The temperature-pressure profiles for models H2-H5 arc 
shown up to extremely low pressures (~10 -8 bars), at which 
point noncquilibrium processes, such as photochemical disso- 
ciation and UV opacity, would be important for determining 
the true thermal profiles. 
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classes of models that fit the five data points op- 
timally. They find that a blackbody temperature 
of —3175 K is needed to fit the Kepler data, sig- 
nificantly higher than the brightness temperature 
inferred by Borucki ct al. (2009); lower brightness 
temperatures are inferred for the IRAC data. Their 
best-fit (non-blackbody) models are (by construc- 
tion) not in local radiative or chemical equilibrium, 
but all of the best-fitting models contain thermal 
inversions, as do our models H2-H5. By relaxing 
the fit-criterion to 1.25-er at each datum and allow- 
ing strongly nonequilibrium chemical abundances, 
they find a noninvcrtcd model with significant CH4 
abundance, but the high temperatures of the atmo- 
sphere favor CO and therefore favor the inverted 
models, consistent with our analysis. 

4. TRES-2 

TrES-2 is a 1.20Mj, l.22Rj planet in nearly cir- 
cular orbit 0.0356 AU from its 1.06 M , 1.00 R Q 
solar-type (GOV) star (O'Donovan et al. 2006). At 
its substellar point, TrES-2 experiences an irradiat- 
ing flux of — 1.1 xlO 9 erg cm~ 2 s _1 , similar to that 
experienced by HD 209458b. 

O'Donovan et al. (2010) report IRAC observa- 
tions of TrES-2, together with nonequilibrium atmo- 
sphere models of the planet generated in the manner 
of Madhusudhan & Seager (2009). O'Donovan et al. 
find that the infrared data can be reasonably well 
fit by a blackbody model, a model with a thermal 
inversion, and a model without a thermal inversion. 
They point out that their model without a ther- 
mal inversion requires a surprisingly low abundance 
of CO, given the temperature of the atmosphere 
(—1500 K). Therefore, they favor the model with 
the inversion. 

Here, we consider a variety of models of TrES- 
2. Motivated by the analysis of O'Donovan et al. 
(2010), we examine models both with equilibrium 
and with nonequilibrium chemistry, all of which are 
in radiative equilibrium. In Table 1, we list five 
TrES-2 models with opacities defined by equilib- 
rium chemistry at solar abundances (T1-T5) and 
five that are completely analogous, except with the 
CO abundance artificially set to (Tln-T5n). In 
the latter group of models, the carbon that would 
have been in CO is instead in CH4, and the excess 
oxygen is instead in H2O. Among both the set of 
models with and without CO, there is a model that 
has no extra absorber and no inversion (Tl and 
Tin, respectively), and four models that do have 
an extra absorber and thermal inversions. After we 
had generated these atmosphere models, Croll et al. 
(2010) and Kipping & Bakos (2010) published Ka- 
band (—2.2 /urn) and Kepler-band observations, re- 
spectively, of TrES-2 's secondary eclipse. The mod- 
els presented herein can, therefore, be thought of 
as predictions, not "postdictions," for the recent 
data. As a result, it was gratifying to see that the 
new data are consistent with our predictions, as de- 



scribed below. 

Figures 2 and 3 portray properties of the at- 
mosphere models with and without CO, respec- 
tively. In these figures, similar to Fig. 1, the top- 
left panel shows the model planet-star flux ratios in 
the optical and the top-right panel shows the same 
in the infrared, while the bottom panel shows the 
temperature-pressure profiles. 

We consider first the equilibrium chemistry mod- 
els in Fig. 2. The non-inverted model (Tl) badly 
fails to reproduce the IRAC data at 4.5 nm and 
at 8.0 /im. The models with inversions (T2-T5) 
have similar temperature-pressure profiles to one 
another. The main difference among these mod- 
els is the temperature in the region of redistribu- 
tion (10~ 2 -10 _1 bars), where the temperature dip 
ranges between -200 K (T2; P n = 0.1) and -700 K 
(T4, T5; P n — 0.3). As a consequence of their simi- 
larity, these four profiles correspond to very similar 
flux ratios at 4.5 ^m, 5.8 /xm, and 8.0 /im. Model 
T5 (k' = 0.3 cm 2 g _1 ) is the best fit at 3.6 /xm, 
but the inverted models all overpredict the 5.8 /Pin- 
point by ~2<t. Models T2-T5 all predict Ka-band 
flux consistent with the Croll et al. (2010) observa- 
tions. All five models predict very low planet-star 
flux ratios in the optical, in contrast to the models 
and the Kepler observation of HAT-P-7b (for which 
F p /F,. = 1.3 x 10~ 4 ; sec §3). Aside from model Tl, 
which is clearly disfavored by the IRAC data, the 
other models all predict planet-star flux ratios of 
<, 2 x 10~ 5 . With no extra absorber, model Tl still 
predicts a Kepler-band flux ratio of only — 3xl0~ 5 . 
Models T2-T5 are all consistent with the Kipping 
& Bakos (2010) observation of (1.1 ± 0.7) x 10~ 5 
(2-er errors). 

The models without CO opacity (Tln-T5n), 
portrayed in Fig. 3, are qualitatively quite similar 
to the equilibrium models in all our diagnostics (op- 
tical flux, infrared flux, and thermal profile). How- 
ever, there are some differences in detail. The deep 
isothermal layers of all five models are hotter by 
—150 K than their equilibrium counterparts, and 
the upper atmosphere of the inverted models (T2n- 
T5n) is cooler by a comparable amount, while the 
upper atmosphere of the non-inverted model (Tin) 
is still warmer than its cqulibrium analog (Tl). The 
magnitude of the temperature dip in the redistribu- 
tion range is also somewhat muted compared with 
the equilibrium models. These slight T-P profile 
differences, in conjunction with the altered opaci- 
ties, result in slightly different optical and infrared 
spectra. In the infrared, the inverted no-CO mod- 
els produce lower flux at 5.8 fim, overpredicting 
the observed data by less than the equilbrium mod- 
els. Among no-CO models, the IRAC data are best 
fit by a model with slightly less extra optical ab- 
sorption (T4n, k' = 0.2 cm 2 g _1 ). These models 
have optical fluxes that are qualitatively similar to 
their equilibrium counterparts, all predicting Ke- 
pler band flux ratios of ^ 5 x 10~ 5 . In particular, 



5 



model Tin predicts slightly higher optical flux than 
does model Tl, but the inverted models T2n-T5n 
predict slightly lower optical fluxes than do models 
T2-T5. The predicted Ks-band and Kepler-band 
fluxes of models T3n-T5n are all consistent with 
observations of Croll et al. (2010) and Kipping & 
Bakos (2010), respectively. 

The basic conclusion that we draw from ex- 
amining models with reduced CO opacity is that 
even such a drastic reduction of CO abundance 
as entirely eliminating it does not cause dramatic 
changes in thermal profiles or spectra. In the 
analysis of O'Donovan et al. (2010), changing the 
[CO]:[CH 4 ]:[C0 2 ] ratio from 10~ 6 : 10~ 6 : (i.e., 
1:1:0) to 10~ 4 : 5 x 10~ 5 : 2 x 10~ 6 (i.e., 50:25:1), to- 
gether with changing the temperature-pressure pro- 
file from a non-inverted to an inverted one, results 
in a modest improvement in the quality of the fit 
(particularly at 8-^.m). We note that, in addition to 
having a surprisingly noncquilibrium [CO]:[CH4] ra- 
tio, the O'Donovan et al. (2010) noninverted model 
has a very low total abundance of carbon. Even 
though their inverted model is substantially sub- 
solar in carbon (by more than an order of magni- 
tude), reducing the carbon abundance by another 
two orders of magnitude is nearly tantamount to re- 
moving carbon from the opacities entirely. In sum, 
our analysis shows that, among models with solar 
abundance of carbon, the presence or absence of CO 
in the database does not make a large difference 
in the thermal profiles or in the emergent spectra, 
and in either case way a hot upper atmosphere and 
thermal inversion are required in the model in or- 
der to come at all close to matching the IRAC data. 
Individual molecular abundances substantially dif- 
ferent from what one would obtain for solar elemen- 
tal abundances might allow for marginally improved 
fits to the data, but might not yet be called for by 
the relatively sparse data available so far. 

The upper left panels of both Figs. 2 and 3 both 
show that models with extra optical absorbers in 
the upper atmosphere have lower planet-star flux 
ratios in the Kepler band than the models without 
an extra absorber. We revisit this point in §5. In 
light of this generic trend, and since among radia- 
tive equilbrium models the IRAC data are better 
fit by inverted models than by non-inverted models, 
we predict that Kepler photometry of TrES-2 will 
indicate low optical flux from this planet, not more 
than one or a few times 10~ 5 of the stellar flux. If 
future Kepler-band observations reveal optical flux 
in excess of this amount, that might be indicative 
of extra optical scattering opacity that was not in- 
cluded in our models. 

5. OPTICAL FLUX VS. OPTICAL ABSORBER 
STRENGTH 

Here, we point out a puzzle: the upper left panel 
of Fig. 1 shows a different trend of optical flux vs. k' 
from the analogous panels of Figs. 2 and 3. Lopez- 



Morales & Seager (2007) noted that inverted mod- 
els, with their hot upper atmospheres, might be ex- 
pected to have higher optical flux than non-inverted 
models. But how do we explain the trend seen for 
TrES-2 (Figs. 2 and 3)? 

If there is an extra absorber in the optical part 
of the spectrum (a K e , in the terminology of Bur- 
rows et al. 2008b, or a k' in the present work), the 
emergent optical flux is affected by two competing 
effects. The absorber makes the planet darker at 
altitude, but can also heat the upper atmosphere. 
The latter effect makes the upper atmosphere of the 
planet more emissive. 

Figure 4 illustrates how the balance of these two 
effects depends on the opacity (n') of the upper- 
atmosphere absorber, in both HAT-P-7b and TrES- 
2. The Kepler bandpass brightness of P n = 0.0 
models of both planets (HAT-P-7b: green; TrES-2: 
blue) is plotted as a function of k', for a series of val- 
ues of k' between and 1.1 cm 2 g _1 . For both plan- 
ets, a small amount of extra absorption results in 
reduced emergent flux in the Kepler band, and, for 
both planets, large values of k' ( > 0.3 cm 2 g _1 ) re- 
sult in increased optical emission, as emission from 
the Wien tail of the hot upper atmosphere becomes 
more prominent in the optical part of the spectrum. 
There is, therefore, a generic non-monotonic char- 
acter to the dependence of optical flux on k' . The 
degree of the non-monotonicity, however, is greater 
for HAT-P-7b models than for TrES-2 models. In 
the former, the optical flux for n' J> 0.3 cm 2 g^ 1 is 
greater than with no extra absorption, while in the 
latter, even with n' = 1.1 cm 2 g _1 the optical flux 
is just over half of what it is with no extra absorp- 
tion. The origin of this difference is related to the 
properties of the irradiation. HAT-P-7 is an excep- 
tionally large (1.84i?0) star that is fairly hot (F6, 
6350 K), whereas TrES-2's star is a more ordinary 
solar-type star (1.003i? Q , GOV, 5850 K). As a re- 
sult, the incident optical irradiation at HAT-P-7b 
is more than 4 times as great as on TrES-2. This 
greatly enhanced incident optical irradiation results 
in a much greater sensitivity to the presence of an 
extra optical absorber. 

6. CONCLUSION 

We have presented atmosphere models of HAT- 
P-7b and TrES-2, two of the three K epler field plan- 
ets that were known prior to the start of the Kepler 
mission. We find that the combination of the IRAC 
and Kepler secondary eclipse data for HAT-P-7b, 
with a Kepler-bandpass secondary eclipse ratio of 
~1.3 x 10~ 4 , appear to require an extremely hot 
upper atmosphere, with an extra optical absorber 
that creates a strong thermal inversion and with 
little day-night redistribution. The IRAC data for 
TrES-2 led us us to expect that TrES-2 has a much 
lower planet-star flux ratio ( ^ 2 x 10~ 5 ) in the Ke- 
pler bandpass than does the HAT-P-7 system, and 
indeed this is what was seen in recently published 
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Kepler data of this object. 

Furthermore, wc find that there is a non- 
monotonic relationship between k' and a planet's 
day-side emergent optical flux. This non- 
monotonicity highlights the need for multiwave- 
length observations in order to better estimate the 
atmospheric structure. 
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APPENDIX 

PARAMETERIZING REDISTRIBUTION IN A ID MODEL 

There have been several approaches used in the literature for treating the redistribution of day-side irradi- 
ation to the nightside in the context of one-dimensional models. Here, we consider several parameterizations 
and the relationships among them. Early work, including that of our group prior to Burrows et al. (2006), 
used only a single parameter (called "/") to describe redistribution. However, there are at least three effects 
that ought to be included in a description of how an atmosphere redistributes heat: (1) some fraction of 
incident energy is transported to the nightside via atmospheric motions (we call this fraction "P n "); (2) 
this redistribution occurs at some depth between the top and the bottom of the atmosphere; and (3) the 
visible face of the planet (the dayside at secondary eclipse phase) has (in general) an anisotropic distribution 
of specific intensity in the direction of Earth. Our more recent work, including this paper, employs our 
attempt, however imperfect, to incorporate these physical effects that ID models from other groups have 
not fully included. In particular, when other modeling efforts have implemented schemes for day-night heat 
transport, they have essentially taken the redistribution to occur before the incoming radiation reaches the 
top of the atmosphere. In contrast, we specify the range of pressures at which the redistribution occurs; 
while our numerical choice might be not be accurate in detail, it is a physically motivated, less ad hoc way 
to parameterize the physics that we know affect emergent spectra. 

One of the best-known redistribution parameters is the afforcmentioned geometrical / factor, which arises 
from an energy balance relation of the following form: 

L p = jRp<rT^ nil , (1) 

or, £ = '(v) 2 - P) 

In eq. (1), L p = (irRp)L*/ (Ana 2 ) is the total stellar power intercepted by the planet (and, therefore, approxi- 
mately equals the emergent radiation from the planet), R p is the planet's radius, L» is the stellar luminosity, 
a is the orbital separation, a is the Stcfan-Boltzmann constant, and T unl { is the temperature of a uniform- 
temperature sphere of radius R p whose flux at Earth near secondary eclipse phase would be the same as 
that of the planet. In eq. (2), F p is the integrated flux from the planet at Earth, 4 and F* is the integrated 
stellar flux. In cqs. (1) and (2), and in the remainder of this section, we have ignored both scattering in the 
planet's atmosphere and its own intrinsic luminosity (from its heat of formation, tidal heating, etc.). The / 
factor has been used by a number of authors in the last decade (Burrows et al. 2003, 2005; Fortney et al. 
2005; Lopez-Morales & Seager 2007; Hansen 2008; and referred to as a in Barman et al. 2005). 5 Perfect 
redistribution (implying that the planet is itself of uniform temperature) corresponds to / = 1/4. Zero 

4 F p = aT^ nii (Rp/d) 2 , where d is the distance from the planet to Earth. 

5 We note that Burrows et al. (2008a) also use a parameter called /, but this / (which is typically set to 2/3) represents the 
direction-cosine of incident irradiation that is used in the planar atmosphere calculation, and is therefore related to distribution 
of temperature on the dayside, instead of to the fraction of energy that is redistributed to the nightside. Some work (e.g. that 
of Barman et al. 2005) integrates the contributions to the total planet-flux at Earth of a series of concentric annuli, from the 
substcllar point to the terminator. Burrows et al. (2008a) show that this integral is extremely well-approximated by taking the 
entire visible hemisphere as being irradiated by the ray at direction-cosine 2/3. 
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redistribution, so that each annulus a given angle away from the substellar point absorbs and reradiates its 
local irradiation isotropically, corresponds to / = 2/3. 6 There is, thus, a "beaming factor" that increases / 
by a factor of 4/3 over the value (1/2) that it would have if the dayside were of uniform temperature instead 
of peaked toward the center of the disk in accordance with the irradiating flux. / is typically implemented 
in ID atmosphere models as simply a uniform reduction of the incident flux at the top of the atmosphere 
(algorithmically, though not physically, equivalent to the effect of albedo). 

Alternatively, we may characterize the redistribution by removing a fraction P n of the incident stellar 
flux on the dayside in a prescribed pressure interval (as described in §2), and inserting the energy at a similar 
(or different) level on the nightside. In such a prescription, P n plausibly ranges between (corresponding 
to no redistribution) and 0.5 (corresponding to the nightside radiating as much power as the dayside, via 
redistributive winds). Burrows et al. (2006), Burrows et al. (2008a), and other work from our group imple- 
ment this procedure. In order to calculate theoretical secondary eclipse fluxes and spectra, one must have a 
model of the three-dimensional distribution of temperature in a planet's atmosphere (e.g., Showman et al. 
2008, 2009; Burrows et al. 2010). In lieu of performing three-dimensional dynamical calculations, our ID 
atmosphere models assume the same beaming factor for day-side emergent radiation, regardless of the total 
day-side radiance; and they assume uniform temperature on the nightside. In particular, on the dayside, 
when P n = 0.5, / = 1/3 (which is 4/3 of the 1/4 that / would equal for a uniform-temperature dayside, 
when half the incident flux has been redistributed to the nightside). More generally, the integrated model 
day-side flux from the planet can be approximately described by cq. (2), taking / to be defined by 

f=\{l-P n ). (3) 

The descriptions in Cowan et al. (2007) of how "Tday" an d "Tnight" depend on P n correctly quantify the 
total emission from the day and night sides of the planet, but they do not include the beaming caused by an 
anisotropic temperature distribution on the dayside that is peaked at the substellar point; therefore, unlike 
what is suggested in Cowan et al. (2007), the integrated day-side flux is not described by taking the dayside 
to be of uniform temperature T<j a y 

Yet a third redistribution parameter, e, is suggested by Cowan & Agol (2010). s ranges between (no 
redistribution) and 1 (full redistribution). It initially appears to be defined similarly to P n , and is described 
(in part) by the following relation: 

where Td ay is the same as T un jf of eq. (1), i.e., a measure of the flux in the direction of Earth. In cq. (4), there 
is a tilde above the expression for day-side temperature to distinguish it from the corresponding expression 
in Cowan et al. (2007), which is a measure of the total day-side emission, as opposed to the flux in the 
direction of Earth. Comparing eqs. (1) and (4) shows that e is a rescaling of /: e = (8 — 12/)/5. Cowan & 
Agol suggest that their e — limit produces a brighter dayside than does the P n = limit, but it does not. 
Instead, the P n = 0.5 limit produces a brighter dayside than does the e = 1 limit (brighter by the beaming 
factor of 4/3). 

Finally, we emphasize that while eqs. (2) and (3) describe how the integrated planet flux depends on 
P n , this parameter does more than simply affect the top-of-atmosphere energy budget. The real benefit 
of P n (which is not captured by the above equations) is that, when nonzero, instead of reducing incident 
irradiation, it removes energy from the dayside (and deposits it on the nightside) at more realistic levels in 
the atmosphere. This process contributes to a slight thermal inversion by cooling the middle atmosphere, 
and, therefore, affects the ratio of 3.6-/im flux to the 4.5-/xm flux in a way that is not reproduced by the 
standard / parameterizations in which incident irradiation is reduced at the top of the atmopshere. 



6 In this case, the specific intensity in the direction of Earth from a annulus at direction-cosine fi away from the substellar 
point (at full-moon phase) is proportional to jj,. Specifically, ![[/,] = [iL* / (2-ko) 2 . 
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TABLE 1 
HAT-P-7B and TrES-2 models 



Model 


^ >n , f i. Chemistry 
(cm 2 g x ) 


Hl a 

H2 

H3 

H4 

H5 


0.10 0.0 equilibrium 
0.05 0.4 equilibrium 
0.10 0.7 equilibrium 
0.10 1.1 equilibrium 
0.00 1.1 equilibrium 


T1 b 

T2 
T3 
T4 
T5 


0.2 0.0 equilibrium 
0.1 0.2 equilibrium 
0.2 0.2 equilibrium 
0.3 0.2 equilibrium 
0.3 0.3 equilibrium 


Tin 
T2n 
T3n 
T4n 
T5n 


0.2 0.0 no CO c 
0.1 0.2 no CO 
0.2 0.2 no CO 
0.3 0.2 no CO 
0.3 0.3 no CO 



a Modcls H1-H5 have the following parameters - orbital semimajor axis (a), planet and stellar surface gravity (g p , g t ), 
planet and stellar radii (R p , R*), and stellar effective temperature T» - appropriate to the HAT-P-7 system (a = 0.0377 AU; 
log 10 g p = 3.38, log 10 g t = 4.08 where g p , are in cgs; R p = 1.363-Rj, where Rj = 7.15 X 10 9 cm is Jupiter's radius; 
i?* = 1.84_Rq; T* = 6350 K). The irradiating spectrum is interpolated from Kurucz models. 

b Models T1-T5 and Tln-T5n have parameters appropriate to the TrES-2 system (a = 0.0356 AU; log 10 g p = 3.30, log 10 g* = 
4.43 where g p , g* are in cgs; R p = 1.224R./; R, = 1.003-Rq; T* = 5850 K). 

c Modcls Tln-T5n have opacities defined by out-of-cquilibrium chemistry, where CO has been artificially removed, with the 
carbon instead in CH4 and the excess oxygen instead in H2O. 
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Fig. 1. — HAT-P-7b model atmospheres. The models (HI— H5) encompass varying degrees of redistribution (parameterized 
as P n ) and upper-atmosphere optical absorption (parameterized as ft', where the values are in cm 2 g — 1 ), and are listed in 
Table 1. Measured data are shown where available. Top left: Optical planet-star flux ratios. The data point (indicated by the 
l-tT error bars) is from Borucki et al. (2009), and the horizontal black line indicates the full width at 10% maximum of the 
Kepler bandpass. The filled colored circles show the integrals of the model flux over the Kepler bandpass. The model that fits 
the data best is H5, which has significant upper-atmosphere absorption (k' = 1.1 cm 2 g" 1 ) and no redistribution. Top Right: 
Infrared planet-star flux ratios. The data points with 1-a error bars are IRAC data from Christiansen et al. (2010). Filled 
colored circles show the integrals of models over the response functions for the IRAC bands. Among models with significant 
upper-atmosphere absorption (models H2-H5), there is little difference in the infrared spectrum. Model H5 is a slightly better 
fit than H2, a slightly worse fit than H3 and H4 at 3.6 /jm, and is a slightly worse fit than models H2-H4 at 4.5 /im. The models 
fit the 5.8-/tm and 8.0-/im IRAC data equally well. Bottom: Temperature-pressure profiles. The models with extra optical 
absorption in the upper atmosphere (models H2-H5) all show large thermal inversions, with upper-atmosphere temperatures 
~1500 K or more hotter than what would be expected in the absence of extra heating (indicated by model HI). Note that our 
models with an extra absorber extend to very low pressures (~10 — 8 bars), where other physics in the thermosphere might be 
important in determining the actual radiative equilibrium thermal profile. 
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Fig. 2. — TrES-2 equilibrium chemistry model atmospheres, analogous to Fig. 1, but for TrES-2 instead of HAT-P-7b. The 
models (Tl— T5) encompass varying degrees of heat redistribution (P n ) and upper-atmosphere optical absorption («;', values in 
cm 2 g — 1 ), and are listed in Table 1. Top left: Optical planet-star flux ratios (log scale). The horizontal black line indicates the 
full width at f 0% maximum of the Kepler bandpass, and is at the level of the data (shown with with 2-cr vertical error bars) 
from Kipping & Bakos (2010). The filled colored circles show the integrals of the model flux over the Kepler bandpass. Models 
with upper atmosphere absorption (T2— T5) tend to have low optical flux (Fp/F* < 2 X 10~ 5 ). Even with k' = cm 2 g — 1 
(model Tl), the optical flux ratio is ~3xl0 -5 . Models T2-T5 are all consistent with the data, to 2 a. Top Right: Infrared 
planet-star flux ratios. The black square data points with 1-cr error bars are IRAC data from O'Donovan et al. (2010), and 
the black cross data point (with l-cr error bars) represents ifs-band data from Croll et al. (2010). Filled colored circles show 
the integrals of models over the response functions for the IRAC bands (at 3.6, 4.5, 5.8, and 8.0 fim) and for ifs-band (at 
2.2 /im). The 4.5-fim, 5.8-/xm, and 8.0-^tm IRAC data clearly favor models with some extra absorption. At 3.6 /im, 4.5 /im, and 
8.0 (im, model T5 is a good fit to the data, but it misses the data by ^2<r at 5.8 fim. Bottom: Temperature-pressure profiles. 
The models with extra optical absorption in the upper atmosphere (models T2-T5) show moderate thermal inversions, with 
upper-atmosphere temperatures ~500 K hotter than what would be expected in the absence of extra heating (indicated by 
model Tl). 
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Fig. 3. — TrES-2 nonequilibrium chemistry model atmospheres, analogous to Fig. 2, but with CO artificially removed from 
the opacities and with CH4 and H2O correspondingly enhanced. Below, we compare with the models shown in Fig. 2. Top left: 
Optical planet-star flux ratios, with Kipping & Bakos (2010) data (including 2-a error bars) shown. As in Fig. 2, models with 
upper atmosphere absorption (T2n— T5n) tend to have low optical flux (F p /F* < 2 X 10 — 5 ). With re' = cm 2 g — 1 (model Tin), 
the optical flux ratio is ~5xl0~ 5 , higher than the corresponding model in Fig. 2, but still low compared with the HAT-P-7b 
data and models. Models T2n-T5n are all consistent with the data, to 2 a. Toil Right: Infrared planet-star flux ratios. Again, 
the IRAC data clearly favor models with some extra absorption. The 5.8-/xm-discrepancy for model T5n is significantly reduced 
compared with the corresponding model with equilibrium chemistry. With CO removed, model T4n is a reasonably good fit 
to the data (within ~ltr at all IRAC bands). Bottom: Temperature-pressure profiles. The overall features of the profiles are 
similar, but the detailed shapes of the profiles are somewhat different, with isothermal regions ~100 K warmer and upper 
atmospheres ~100 K cooler than the equilibrium chemistry models. 
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Non-monotonicity in Kepler-band- 
integrated flux vs. opacity (k) 
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Fig. 4. — Kepler bandpass brightness ratios vs. upper atmosphere absorption (k'). Models of HAT-P-7b (thick light green) 
and TrES-2 (thick blue) are shown, in addition to the Kepler data for HAT-P-7b (dashed red) with 1-c errors (dotted dark 
green). These models have no redistribution to the nightside (P n = 0). The relationship between optical flux and n' is non- 
monotonic for both planets, but more so for HAT-P-7b. In particular, for HAT-P-7b, large enough values of k' (values above 
~0.3 cm 2 g —1 ) lead to enough flux in the region of the Wien tail that the optical flux is greater than it would be with no 
extra absorber. However, for TrES-2, despite the non-monotonicity, the planet is still brightest in the optical without any extra 
absorber. 



